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^ulkHeatFluxParameterizatio^ 

•  Total  heat  flux  is  available  from  archived  products 

o  SST  drifts  if  total  flux  alone  is  used  to  force  an  OGCM 

•  HYCOM  uses  model  SST  and  a  bulk  heat  flux  parameterization 

•  Feedback  between  SST  and  heat  flux  to  prevent  SST  drift 

•  No  need  for  explicit  relaxation  to  SST 

•  MICOM  uses  constant  exchange  coefficients 

•  HYCOM  has 

o  several  options  for  exchange  coefficients 
o  a  blackbody  longwave  correction 


Subsurface  Heating  Parameterization 


Net  heat  flux  at  a  given  depth,  z: 

Q{z )  =  Q(0)  +  [Qsw(O)  —  Qsvr{z)], 

Net  heat  flux  absorbed  at  the  sea  surface,  z  =  0: 

Q(0)  =  Qlw  +  Ql  +  Qs: 

o  Qlw-  net  longwave  radiation  at  the  sea  surface, 
o  Qsw-  net  shortwave  radiation  at  the  sea  surface, 
o  Ql\  latent  heat  flux, 
o  Qs :  sensible  heat  flux. 

NOTE: 

o  HYCOM’s  “surface”  heat  flux  is  not  Q(0),  but 
o  rather  the  near  surface  flux  absorbed  in  layer  1 
o  e.g.,  Q(  1)  when  the  top  model  layer  is  1  m  thick. 


(1)  Longwave  Radiation 


•  Input  from  archived  products  (e.g.,  ECMWF,  NCEP,  etc) 

•  A  correction  is  needed.  Why  ? 

o  They  use  their  model  SST 
o  Different  from  HYCOM  SST 

•  HYCOM  uses  a  blackbody  correction  (Kara  et  al  2004a) 

Qun(Ts)  =  Qhw{Tsa )  —  (4.506  —  0.0554  TS)(TS  —  Tsa). 

o  Ts:  HYCOM  SST 
o  Tsa :  Atmospheric  model  SST 

•  The  effects  of  clouds  are  independent  of  SST 


-2  o 


A  linear  approximation  to  the  blackbody  radiation 


•  Cubic  formulation  (Josey  et  al.  2003) 

•  Linear  approximation  (Kara  et  al.  2004a) 


(2)  Shortwave  Radiation 


Previous  parameterizations  in  HYCOM 
o  Jerlov  water  types  (Halliwell  2004) 
o  2-band  scheme  (Kara  et  al.  2004b) 

o  Turbidity-dependent  split:  red  and  blue  spectrums. 
o  Based  on  SeaWiFS  /cpar  climatology  (2004c). 
o  Attenuation  coefficient,  &par:  depth-independent 

New  parameterization  in  HYCOM  (in  progress) 
o  Fixed  frequency  ranges: 

o  visible  spectrum  (350-700  nm),  also  called  PAR 
o  infrared  spectrum  (700-2400  nm) 
o  Will  use  absorption  and  backscattering  coefficients 
o  A:par  depends  on  depth  and  solar  angle 


The  shortwave  radiation  at  a  given  depth  (z)  is  split  into  two  parts: 


Qsw(^)  =  Qpar(^)  +  Qm{z),  (3) 

Qpar(^)  =  Qpar(O)  exp(— z  kpAR),  (4) 

Qm(z)  =  Qm(0)exp(-z  fcm),  (5) 


Qir(z)/Q[R(0)  ==>  700  to  2500  nm 
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<2PAR(z)/(9pAR(0)  :  350  to  700  nm 
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(3)  Latent  Heat  FluxTj 


Bulk  formulation:  QL  =  pa  CL  L  Va  (qa  -  qs) 

Previous  exchange  coefficient  (Cl)  in  HYCOM  (Kara  et  al.  2002) 
o  based  the  COARE  (v2.6)  algorithm  (Fairall  et  al.  1996) 
o  excluded  Va  <  4  m  s-1,  Va  >  20  m  s-1 
o  Cl  was  dependent  on  (Ta-Ts)  and  Va 

New  Cl  parameterization  in  HYCOM  (Kara  et  al.  2004d) 
o  based  on  the  COARE  (v3.0)  algorithm  (Fairall  et  al.  2003) 
o  includes  Va  from  1  to  40  m  s_1 
o  Cl  is  dependent  on  (Ta-Ts),  Va ,  and  RH  as  well 

NOTE:  Calculate  Ql  using  HYCOM  SST  at  each  time  step 


Previous  exchange  coefficients  for  the  latent  heat  flux  (RH=100%) 
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New  exchange  coefficients  for  the  latent  heat  flux  (for  varying  RH) 
Is  including  relative  humidity  in  Cl  important? 
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Air-sea  temperature  difference  (  C) 


What  could  be  the  typical  error  in  latent  heat  flux  without  RH? 
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Latent  heat  flux  coefficient  for  RH=100%  (103  CL) 


Latent  heat  flux  (W  m  2) 
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Ql  =  pa  Cl  L  Va  ( qa  ~  qs)  ,  where  Ta  -  Ts  =  2°C,  qa  -  qs  =  3  g  kg  1 


Summary 


•  Global  HYCOM  simulation  with 

o  the  RH-dependent  exchange  coefficients 
o  the  depth-dependent  shortwave  radiation 

•  Shortwave  radiation  attenuation: 

o  need  two  satellite-based  input  fields  for  HYCOM 
o  (1)  absorption  coefficient 
o  (2)  backscattering  coefficient 
o  form  a  climatology  (2001-2003)  using  MODIS 
MODIS:  Moderate  Resolution  Imaging  Spectroradiometer 


